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Abstract 

Elastic singularities such as crack tips, when in motion through a medium 
that is itself vibrating, are subject to forces orthogonal to the direction of 
motion and thus impossible to determine by energy considerations alone. 
This fact is used to propose a universal scenario, in which three dimension- 
ality is essential, for the dynamic instability of fast cracks in thin brittle 
materials. 

The dynamics of crack tips is of particular relevance to the fracture of materials [|TJ 
and to earthquake rupture ||. In addition, recent improvements in instrumenta- 
tion techniques have revealed unstable crack dynamic behaviour reminiscent of 
what is observed in other physical systems driven far from equilibrium ||, f|, [5|. |(|. 
These experiments have taken a close quantitative look at the crack-tip velocity 
fluctuations, ultrasound emission and surface heterogeneities in glass and plexi- 
glass associated with long standing puzzles in the tensile cracking of very brittle 
solids [0]. 

The theoretical effort directed at an understanding of these phenomena has 
been largely based on two dimensional modelling. It includes analysis based 
on a continuum approach that considers a cohesive region at the crack tipf| 
and numerical computations using lattice as well as molecular dynamics ||10|| 
techniques. However, ultrasound [|, ||, surface topography || || and velocity 
fluctuation p, [5|] measurements clearly involve length scales comparable or smaller 
than sample thickness, strongly suggesting that three dimensionality cannot be 
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ignored. The possible relevance of results obtained in a study of the three dimen- 
sional propagation of a crack within an heterogeneous model elastic solid have 
been recently pointed out l|il"H . 

The purpose of this paper is to point out that, when the motion of an elastic 
singularity deviates from a straight line, energy balance considerations are in- 
sufficient to determine the singularitiy's motion since (in two dimensions) they 
provide only one equation to determine the dynamics of two degrees of freedom. 
Thus, in addition, the equation of momentum balance must be considered. In 



so doing there results a new elastic force that does no work[|12] and that conse- 
quently cannot be obtained on the basis of energy considerations alone. For a 
crack moving in a thin elastic plate, consideration of three dimensionality makes 
this force operative and provides a scenario, some of it qualitative, some of it 
quantitative, that appears to explain many of the observed experimental facts. 

The equations of dynamic elasticity in a medium of dimensionality d and 
density p can be obtained as extrema of the action functional 
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where is the stress, Sij the strain 



1, . . . , d and u(x, t) the particle 
deviation at time t from an equilibrium position x. When stress and strain are 
linearly related, = CjjjySjy, the well known linear equations can be written in 



conservation-law form 13, 141: 
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The components of the energy momentum tensor T have the following interpre- 
tation: T 00 is the energy density, T - 7 is the energy flux, T j0 is the momentum 
density and T J * is the momentum flux. 

If Eqns. (0) are integrated over a volume V that is bounded by a surface S 
that is possibly moving with velocity V one obtains energy conservation 
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J d d xT 00 = J dSjT 0j - J dSjV^T 00 (8) 
and (linear) momentum conservation 
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If Eqn. (|8]) is applied to the case of a crack advancing in a straight line, we get 
the usual equation of motion for the crack tipQ. There is energy flow into the 
tip, which is equated to the surface energy that is created as the crack advances. 
The volume V is the volume of the whole sample, and the surface S is composed 
of a small (i.e. smaller than any relevant length such as sample thickness or 
radiation wavelength) cylinder surrounding the crack tip, the lips of the crack, 
and the boundaries of the sample. The contribution of the small cylinder around 
the tip is finite because stresses behave like 1/y/r, with r distance to the tip, 
near the crack. The crack lips do not contribute because they are at rest and, by 
definition, they are traction free. Remote surfaces do not contribute as long as 
radiation reaction is neglected. 

It is a simple matter to check that, for a crack advancing along the x 1 direction 
at velocity V the following is true: 
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E = [ d d xT 00 (11) 
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V 

so that the momentum equation and the energy equation provide the same in- 
formation. Now, if the crack tip motion is not constrained to a straight line, the 
momentum equations provide the additional information needed to determine the 
motion. In two dimensions, one equation would not suffice and less so in three. 

Consider now the normal modes of vibration of an elastic plate in three 
dimensions ||15|j. For a thin plate of thickness h there are, to leading order, one 
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dilatational mode with frequencies nna/h and two shear modes with frequencies 
mr/3/h, where a is the speed of sound and (3 the speed of shear waves, both in 
bulk, and n = 1,2,3,.... (Figure 1). When a crack propagates along the x l 
direction, say, because of static loading along a; 3 , the modes involving motions 
parallel to the crack faces will be unaffected. There is, however, one shear mode 
that involves motion along the x 3 direction that will be incompatible with the 
traction-free boundary condition that defines a crack. Consequently this mode 
will be modified by the presence of the crack and stress concentration, in addition 
to the one already present due to the static loading, will result. The question 
naturally arises as to how this new stress concentration at the crack tip will affect 
its dynamics. Of course, these vibrating modes do not exist in two dimensions. 

In order to give a quantitative assesment of this effect we consider the crack 
surface near the tip as a continuous distribution of infinitesimally small disloca- 
tions. These dislocations are defects in a continuum, and they are not related 
to a (possibly non-existent) underlying crystal structure. That this makes sense 
is well-known in the static case [|L6|] and has also been used in dynamics, to show 
that a crack may, under certain circumstances, acquire inertia[]l?|| . The basic 
idea is that the infinitesimal dislocations arrange themselves so as to negate the 
externally applied stress and enforce the no traction boundary condition at the 
crack surface. The exact functional behaviour of the crack opening displacement 
that will result in the case at hand is not needed at the current stage of develop- 
ment of the theory, although it presents itself as an attractive challenge for future 
research. We only need to know that some such distribution exists. 

Consider then an edge dislocation of (infinitesimal) Burgers vector b moving 
with velocity V in an elastic medium whose (small) particle displacements away 
from an equilibrium position are U(x,t). The relative orientation of b and V 
is immaterial. The time dependence of U does not include a possible motion 
of the body as a whole, since that would not be a small deviation away from 
an equilibrium position. The dislocation generates around it a well known [15] 
particle displacement u(x, t) in which strain and particle velocity are related, at 
short distances, by (a, c = 1, 3) 

The point is to compute the force 



given by (||) and fll2|) , for which the relevant integrals have to be evaluated around 
a small circle enclosing the dislocation. Assuming the medium to be linear down 
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to this small circle, the particle displacements will be the sum 

u + U (15) 

of dislocation generated plus externally imposed ones. Since the components of 
the energy momentum tensor T are bilinear functions of derivatives of displace- 
ments, they will decompose naturally in three terms: one involving only u, one 
involving only U, and one involving both u and U which is the only one of interest: 

I '(9?;* BIP BIP 8v i \ 



and similarly for T lJ . Substitution of the field u for an edge dislocation |L8[ into 
©, (|I2D and (Jig) using QT6|) gives the result 

F ] =F PK + F^ (17) 

where 

— biY^ijejk m t m (18) 

with t the unit vector along the dislocation line and Ejj = Cijk m dUk/dx m the 
external stress, is the usual Peach-Koehler force. In addition there is a force ||19|| 



F v = pejkmVkimbiUi (19) 

that acts only if the dislocation is moving (V ^ 0), if the elastic medium itself 

is also independently moving (U ^ 0), and is perpendicular to the direction of 
motion and consequently does no work. For a Burgers vector corresponding to 
a mode-I crack, the static loading gives a Peach-Koehler force along the cleav- 
age plane. When summed over the many infinitesimal dislocations that make 
up the crack tip, it gives the usual work-performing force that is considered in 
the standard treatement of crack-tip dynamics . It does not have a component 
perpendicular to the motion. The vibrational modes of the plate that have been 
previously considered give a vanishing Peach-Koehler force. The force Fy, on 
the other hand, is nonvanishing only for the shear mode that was already argued 
to be concentrated at the tip in order to satisfy the proper boundary conditions 
along the crack faces. The static loading of course does not contribute to it. 

As it was mentioned at the beginning of this communication, recent experi- 
ments H f|, |5], Q indicate that cracks in glass and plexiglass suffer an instability 
at a critical velocity, beyond which there are strong sound emissions, strong ve- 
locity fluctuations, and strong surface heterogeneities. The limiting velocity is 
well below the theoretically expected Rayleigh wave velocity. We are now in a 
position to propose the following scenario to explain those observations: 



5 



1. The crack starts propagating, and right away it starts emitting elastic (both 
longitudinal and shear) radiation. The pull tabs possibly also radiate at the 
beginning. The crack accelerates. 

2. Normal modes (in thickness) of vibration of the plate are established, and 
they are responsible for a new oscillating force that does no work on the 
crack and is proportional to the tip's velocity. The fundamental mode will 
dominate but at sufficiently high crack velocity the forcing may be enough 
to also excite higher modes. To leading order they are the modes of an 
infinitely extended slab. Next order corrections should take into account 
the finite size and shape of the sample. 

3. The crack responds by oscillating to the tune of this new periodic (possibly, 
multiply periodic) forcing and, hence, emitting radiation (both shear and 
longitudinal) primarily at the same frequencies it is being forced. This 
radiation reinforces the normal modes that, in turn, feed back on the tip. 
This is the instability. 

4. In a non-amorphous material such as a crystal, the surface energy may be 
higher for out of plane crack motions and thus this mechanism may be non 
operative. 

5. Terminal speed (i.e. average speed along line of propagation) is determined 
by energy spent by external stress. Since the crack path is effectively in- 
creased, the average speed of advance is decreased. 

The above scenario is universal in the sense that it applies to any brittle material 
that is described by linear elasticity. At this stage, there is no calculation of a 
transition velocity, but a reasoning that indicates why there should be an insta- 
bility as the crack velocity increases, There is, however, a clear consequence that 
sound radiation, velocity oscillations and surface roughness are different manifes- 
tations of the same phenomenon, whose spectra should be peaked at the normal 
modes of the elastic plate. These consequences should be amenable to experi- 
mental verification, for example by studying the effect of sample thickness and of 
sample ensonification. 

To conclude: it has been established that elastic singularities are subject to 
velocity dependent forces that do no work. On this basis, a scenario has been pro- 
posed to explain recent observations of dynamic crack instabilities. The proposed 
scenario includes the three dimensionality of the problem as an essential ingredi- 
ent. It provides a qualitative physical origin for the instability, and a quantitative 
prediction for the frequencies that should be notorious in the resulting spectra of 



6 



sound, surface and velocity fluctuations, namely, the eigenfreuencies of the thin 
plate that is being fractured. 
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Caption for figure 

Figure 1 A thin elastic plate of thickness h has, to leading approximation, three 
normal modes of elastic oscillation: one (a) is dilatational, with eigenfrequencies 
nna/h and two (b-c) are shear, with eigenfrequencies nirj3/h, where a is the 
speed of sound, (3 is the speed of shear waves, and n = 1,2,... When a crack 
propagates along the ^direction, for instance because of static loading along 
a; 3 , the shear mode indicated in (c) is incompatible with traction-free boundary 
conditions along the crack faces, leading to a new stress concentration at the tip 
that will modify the crack dynamics. 
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